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ABSTRACT 


An  array  of  three  borehole  tiltmeters  near  Quebec  City  in  Eastern 
Canada  is  designed  to  study  the  tidal  and  secular  response  of  the  crust 
in  the  Charlevoix  seismic  zone.  The  objectives  of  this  first  study  of 
data  from  two  boreholes  of  the  array  re  to  investigate  the  spatial 
coherency  of  the  tidal  observations  and  determine  whether  there  are 
time  variations  in  the  tidal  amplitudes  and  phases;  and  to  explore  the 
main  features  of  the  secular  tilt  signal.  The  tidal  analysis  was  done 
using  a  modified  version  of  the  HYODN  harmonic  analysis  program  of 
Schuller  [1977] ,  with  which  the  time  varying  tidal  amplitudes  and 
phases  were  determined  by  the  sequential  analysis  of  overlapping  2 
monthly  subsets  of  the  data.  The  admittance  variations  observed  for 
the  major  constituents  (M2,  N2,  S2K2r  01,  PlSlKi)  show  remarkable 
agreement  between  boreholes,  suggesting  that  they  are  regional  in 
origin.  Differences  of  up  to  20%  in  anplitude  and  5  degrees  in  phase 
were  found  between  the  mean  M2  results  determined  frcm  boreholes  1  and 
2,  indicating  small-scale  distortion  of  the  local  tilt  field  by  lateral 
inhomogeneities.  The  secular  tilt  frcm  both  boreholes  correlates 
strongly  with  transient  and  seasonal  water  table  fluctuations,  suggest¬ 
ing  the  dominant  influence  of  pore  pressure  effects  on  the  non-tidal 

* 

tilt.  A  preliminary  estimate  of  the  detectability  of  long  term  region¬ 
al  trends  in  tilt  is  0.4  ^irad/yr. 


1.  INTRODUCTION. 


An  array  of  three  borehole  tiltmeters  has  been  established  in  the 
Charlevoix  region  of  Quebec  (figure  1) ,  an  area  considered  to  be  among 
the  most  seismically  active  in  eastern  North  America.  The  tilt  observ¬ 
ations  are  being  made  at  the  Charlevoix  Geophysical  Observatory  of  the 
Canadian  Department  of  Energy,  Mines  and  Resources  and  are  aimed  at 
investigating  the  behaviour  of  crustal  rocks  in  an  earthquake  zone. 
Specifically,  the  program  of  measurements  is  designed  to  saitple  the 
earth  tide  response  as  a  possible  indicator  of  changing  crustal  condit¬ 
ions  and,  by  monitoring  secular  and  transient  tilts,  to  directly  detect 
regional  crustal  deformations  associated  with  processes  occurring  at 
depth. 

This  first  analysis  of  the  borehole  tilt  data  from  the  Charlevoix 
Observatory  provides  an  explanation  of  the  first  order  properties  of 
the  secular  tilt,  establishes  a  methodology  for  determining  the  real 
variations  in  the  tidal  admittance,  and  applies  this  approach  to  the 
Charlevoix  data.  On  the  basis  of  our  results  we  then  assess  the  potent 
ial  of  using  the  time  varying  tidal  admittance  at  Charlevoix  as  an 
indicator  of  changes  in  crustal  response;  and  of  using  the  secular 
component  to  detect  long  term  tilts  of  tectonic  origin. 

The  first  study  of  tidal  response  in  relation  to  earthquake  activ¬ 
ity  was  reported  by  Nishimura  [1950]  who  recorded  fluctuations  of  up  to 
75%  in  monthly  amplitude  determinations  of  the  M2  constituent  using 
quartz  pendulums  installed  165m  below  ground  in  Makimine  copper  mine, 
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Japan.  These  variations  were  attributed  to  elastic  parameter  changes 
in  the  local  crust  modifying  the  tidal  response.  It  was  suggested  that 
the  longer  term  changes  in  the  tidal  amplitudes  were  correlated  with 
fluctuations  in  seismic  activity.  Other  studies  have  indicated  less 
dramatic  effects,  the  measured  tidal  response  variations  generally  not 
exceeding  15%.  Kato  [1979]  measured  2-4%  amplitude  fluctuations  in  M2 
from  observations  recorded  by  a  water  tube  tiltmeter  at  Kamitakara 
during  September,  1977  to  October,  1978.  Monthly  M2  strain  amplitudes 
at  the  same  location  fran  1968  to  1972  [Mikumo  et  al,  1978]  showed  a 
steady  rise  of  15%  during  the  eleven  months  prior  to  the  1969  Gifu 
(M=6.6)  earthquake.  More  recent  results  fran  the  Makimine  station 
[Tanaka,  1976]  indicated  amplitude  variations  not  exceeding  a  range  of 
20%;  and  it  has  been  suggested  that  the  extraordinary  amplitude  changes 
seen  by  Nishimura  [1950]  included  instrumental  effects.  Agnew  [1979] 
has  studied  the  time  variation  of  M2  in  strain  data  recorded  on  the 
800m  laser  strainmeter  at  Pinon  Flat,  California.  Fran  1973  to  1978 
fluctuations  were  confined  within  2%  of  the  mean  amplitude.  Given  the 
high  signal  to  noise  ratio  (45dB)  of  the  estimates  and  the  fact  that 
Pinon  Flat  is  locally  an  aseismic  region,  this  result  establishes  an 
experimental  baseline  upon  which  to  evaluate  the  significance  of  time 
variations  in  admittance  in  active  seismic  zones. 

Beaumont  and  Berger  [1974]  suggested  that  if  Vp/vs  seismic  veloc¬ 
ity  anomalies  are  indicative  of  changes  in  crustal  elastic  properties, 
then  they  should  be  accompanied  by  tidal  response  anomalies.  For  tilt, 
this  would  be  caused  by  a  laterally  varying  effect  where,  for  example, 


a  15%  reduction  in  Vp  within  a  crustal  inclusion  would  be  accompanied 
by  a  maximum  30-40%  change  in  amplitude  at  the  edge  of  the  inclusion, 
decreasing  rapidly  with  distance  away  from  it.  Tanaka  [1976]  has  used 
a  similar  approach  to  look  at  the  modification  of  the  marine  load  tide 
response  due  to  a  dilatant  zone  near  the  coast.  Beaumont  [1978a] 
expanded  on  the  work  of  Beaumont  and  Berger  [1974]  by  considering  the 
effect  on  the  tidal  response  of  the  superposition  of  tidal  and  increas¬ 
ing  tectonic  stress  regimes  in  the  crust.  Appealing  to  laboratory 
results  which  showed  that  hysteresis  occurs  in  intact  rock  samples  as 
they  approach  failure,  he  outlined  ways  in  which  linear  and  non-linear 
anomalies  in  the  tidal  response  may  be  generated  depending  on  the  rate 
of  change  of  stress  in  the  crust.  Agnew  [1981]  extended  the  theoret¬ 
ical  analysis  further  considering  different  types  of  non-linearities 
and  searched  for  their  signature  in  strain  records  at  Pinon  Flat 
without  success. 

There  is  a  catalogue  of  attempts  to  measure  tectonic  tilts  and 
strains  in  the  non-tidal  part  of  the  spectrum  (see,  for  example,  Wyatt 
et  al  [1982];  McConnell  and  Lewkowicz  [1978];  Herbst  [1976];  Mortensen  and 
Johnston  [1975]).  Few  have  been  demonstrably  successful,  for  the  most 
part  because  such  measurements  are  carried  out  near  the  surface  where 
they  are  severely  polluted  by  other,  inseparable  signals  usually  of 
meteorological  origin.  It  is  clear  that  any  hope  of  success  lies  in  an 
approach  which  either  ensures  the  isolation  of  the  measurement  frcm 
these  sources  of  noise  (by  deep  enough  burial)  or  one  which  entails  a 
spatial  average  of  the  near  surface  signal.  This  second  approach  may 


involve  a  network  of  instruments  and  subsequent  statistical  elimination 
of  noise  [Morten9en  and  Johrston,  1975],  or  the  use  of  long  baseline 
instruments  [Wyatt  et  al,  IS  82].  These  methods  assune  that  the  spatial 
scales  of  the  meteorologies]  and  other  sources  of  noise  are  less  than 
that  of  true  tectonic  tilt  cncmalies.  Wyatt  et  al  [1982]  in  an  inport- 
ant  comparative  experiment  et  Pinon  Flat,  demonstrated  the  relative 
potential  for  secular  tilt  Etudies  of  measurements  performed  over  a 
long  baseline  (535m),  at  mocerate  depth  (26m),  and  in  shallow  (5m) 
borehole  installations.  Their  results,  for  data  spanning  a  5  month 
period,  indicate  successive  order  of  magnitude  i)  creases  in  the  drift 
rates  from  the  long  baseline  tiltmeter  (0.07  pr ad/month)  to  the  26m 
deep  borehole  tiltmeter  (0.4  prad/month)  to  the  shallow  installations 
(up  to  7  prad/month) .  However,  the  most  encouraging  results  to  date  are 
those  of  Cabaniss  [1978]  who  achieved  a  net  secular  drift  of  0.3  prad 
in  a  120m  deep  borehole  at  Eedford,  Massachussetts  over  a  three  year 
period. 

The  Charlevoix  region  has  a  history  of  large  earthquakes  occurring 
roughly  every  60  to  90  years.  The  last  major  earthquake  was  a  magnit¬ 
ude  7  which  struck  in  1925.  Seismicity  surveys  in  1968,  1970  and  1974 
[Milne  et  al,  1970;  Leblanc  et  al,  1973;  Leblanc  and  Buchbinder,  1977] 
show  that  the  earthquake  zone  covers  a  70  by  40  Km  region  centred  on 
the  St.  Lawrence  River  150  Km  northeast  of  Quebec  City.  Almost  all  of 
the  earthquakes  were  shown  to  have  occurred  within  the  Precambrian 
rocks  of  the  Grenville  Province,  at  an  average  depth  of  11  Km.  The 
nodal  plane  solutions  are  contpatible  with  pure  thrusting,  the  planes 
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dipping  on  average  40  degrees  to  the  west  and  50  degrees  to  the  east 
[Leblanc  and  Buchbinder,  1977]  with  the  direction  of  the  pressure  axis 
not  in  disagrement  with  what  Sbar  and  Sykes  [1973]  found  in  other 
parts  of  North  America.  Overlying  the  Grenville  basement  rocks  are  the 
Ordovician  sediments  deposited  during  the  Taconic  orogeny.  During  the 
1974  survey  no  events  were  detected  above  the  contact  between  the 
basement  and  sediments,  which  outcrops  at  Logan's  Line  along  the  north 
shore  of  the  river.  It  appears,  on  the  basis  of  this  measured  distrib¬ 
ution,  that  the  length  of  the  seismic  zone  parallel  to  the  river  is 
related  to  the  350Ma  old  Charlevoix  inpact  crater  [Robertson,  1968; 

1975]  (see  figure  1) .  The  conjunction  of  Logan's  Line  and  the  crater 
could  influence  the  seismicity  distribution  [Leblanc  et  al,  1973; 

Leblanc  and  Buchbinder,  1977] .  Presumably,  the  structures  either  repr¬ 
esent  a  zone  of  weakness,  remaining  from  crustal  damage  caused  during 
the  impact,  in  which  tectonic  stress  is  released  by  earthquakes,  or  the 
structures  focus  the  tectonic  stress  in  sane  way  causing  a  higher  rate 
of  seismicity  in  an  otherwise  normal  crust. 

Seismic  crustal  studies  [Lyons  et  al,  1980]  indicate  that  the 
Ordovician  overthrusts  dip  beneath  the  river  southeast  of  Logan's  Line 
at  an  average  of  20  degrees.  Their  results  also  indicate  that  the  upper 
crustal  seismic  P-velocity  in  a  belt  extending  along  the  north  shore  of 
the  river  is  0.2-0. 3  KiVsec  lcwer  than  the  average  6. 4-6. 5  Km/ sec  for 
the  Grenville  Province.  They  suggest  that  this  is  related  to  structural 
deformation  caused  by  the  Charlevoix  meteorite  inpact.  The  near  surface 
P-velocity  within  the  crater  is  6.08  Km/ sec  suggestive  of  a  zone  weak- 
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ened  by  fracturing  during  inpact. 

Since  1972  various  kinds  of  geophysical  measurements  have  been 
made  by  the  Earth  Physics  Branch  of  the  Department  of  Energy,  Mines  and 
Resources  to  supplement  the  existing  horizontal  and  vertical  control 
provided  by  the  Geological  Survey  of  Canada.  Thus,  at  the  present  time 
the  geophysical  measurements  in  the  area  [Buchbinder  et  al,  in  press] 
comprise:  1)  first  order  levelling  along  the  north  and  south  shores  of 
the  St.  Lawrence  River;  2)  a  precise  gravity  network  consisting  of  15 
stations  and  twenty-two  connections;  3)  a  triangulation  network  consis¬ 
ting  of  twelve  points  and  twenty-six  connections  spanning  the  St. 
Lawrence  River;  4)  two  permanent  tide  gauges;  5)  a  permanent  seismic 
observatory  at  La  Pocatiere  on  the  south  shore;  6)  an  array  of  four 
magnetotelluric  stations;  and  7)  tilt,  strain,  seismic  and  magnetotell¬ 
ur  ic  installations  located  at  the  Charlevoix  Observatory  and  operating 
in  conjunction  with  borehole  tiltmeter  programme  to  be  discussed  in 
this  study. 

In  the  following  sections,  we  describe  briefly  the  experimental 
set-up  and  discuss  in  detail  the  methods  and  procedures  used  in  analys¬ 
ing  the  tilt  data.  Particular  attention  is  paid  to  the  elimination  of 
apparent  time  variations  in  the  admittance  which  are  artifacts  of  the 
analysis.  It  is  shown  that  while  the  mean  admittances  for  the  M2  and  Oi 
tidal  constituents  are  in  poor  agreement  between  boreholes,  there  is 
good  agreenent  between  the  boreholes  for  the  identified  real  time 
changes  in  admittance.  Quasi-per iodic  variations  in  the  M2  constit¬ 
uent,  indicating  an  approximately  six  monthly  cycle,  reach  a  level  of 
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5%  in  amplitude  and  3-4  degrees  in  phase.  Non-tidal  tilts  in  the  obser¬ 
ved  data  are  at  the  level  of  1-2  >irad  and  are  strongly  correlated  with 
groundwater  level.  Incoherent  residual  drifts  appear  to  be  present  at 
a  minimum  level  of  0.4  yarad/yr . 


2.  TOE  BOREHOLE  TILTMETER  ARRAY 

The  tiltmeter  borehole  array  consists  of  three  30cm  diameter 
holes,  two  of  which  are  46m  deep  (boreholes  1  and  2)  and  the  third 
(borehole  3)  which  is  110m  deep  [Peters,  1983] .  Each  is  cased  with  a 
20cm  diameter  mild  steel  casing  terminated  by  a  6m  long  stainless  steel 
pod  of  the  same  diameter  which  houses  the  tiltmeter.  All  of  the  holes 
are  dry,  but  only  the  shorter  two  are  straight.  Figure  2  is  a  diagram 
of  the  experimental  configuration,  showing  a  plan  viw  of  the  observat¬ 
ory,  and  the  tiltmeter  installation  in  section.  Adjacent  to  each  of 
the  boreholes  is  a  well  (referred  to  as  A,  B  or  C)  in  which  manual  and 
continuous  measurements  of  water  level  are  made.  A  vault  containing 
three  ANAC  mercury  level  tiltmeters  is  situated  near  the  centre  of  the 
borehole  array  [Peters  et  al,  1983]. 

■in.  ehnle  tiltmeters  are  Bodensewerk  Gbp  10  (formally  Askania) 
vertical  pi  ulums  [Flach  and  Rosenbach,  1971]  chosen  because  of  their 
proven  per  urmance  [Zschau,  1976;  Edge  et  al,  1979;1983;  Flach  et  al, 
1975]  and,  not  least,  because  they  are  the  only  high  quality  tiltmeters 
commercially  available  and  supported.  The  most  important  feature  of  the 
instrument  is  the  accurate  in  situ  calibration  absent  in  most  other 
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ions  between  105  and  107,  operating  at  different  times  in  borehole  1, 
exceeds  the  range  of  variability  so  far  determined.  The  averaged  ANAC 
results  do,  however,  plot  between  the  borehole  1  observations.  Tilt- 
meter  106  in  borehole  2  lags  the  other  measurements  by  4-5  degrees  in 
the  south  direction  and  leads  by  nearly  7  degrees  in  the  east  direct¬ 
ion.  Furthermore,  there  is  a  20%  difference  between  the  106  and  107 
amplitudes  in  the  east.  The  0^  results  are  in  fairly  good  agreement.  A 
detailed  examination  of  the  coherency  of  a  number  of  constituents 
[Peters  and  Kumpel,  in  prep]  has  confirmed  that  the  large  discrepancies 
between  the  two  boreholes  for  seme  of  the  constituents  is  a  real  effect 
rather  than  instrumental. 

While  there  is  general  agreement  between  theory  and  the  observed 
Op  admittances,  there  are  systematic  differences  between  all  of  the  M2 
measurenents  and  theory.  In  both  directions,  observations  lag  the  model 
results.  In  the  south,  where  drying  is  iirportant,  the  observed  amplit¬ 
udes  lie  within  the  range  of  the  total  drying  (3)  and  non-drying  (4) 
model  vectors  and  may  constrain  the  model  to  the  partial  drying  case 
which  is  midway  between  the  two  extremes  and  is  physically  the  most 
realistic.  Nevertheless,  the  layered  earth  model  used  to  compute  the 
tilt  Green's  functions  does  not  take  into  account  the  regional  struct¬ 
ure  dominated  by  Logan's  contact  (section  1)  nor  are  corrections  made 
for  the  effects  of  local  topography  and  geological  structure. 

5.4  Time  variations  in  the  tidal  constituents 


One  aspect  of  this  study  is  to  identify  "true"  variability,  or 
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higher  frequencies  in  Berger  and  Levine  [1974]  .  Although  the  tilt 
spectra  do  not  extend  to  high  enough  frequencies  to  show  a  conclusive 
trend,  there  is  the  suggestion  of  a  similar  increase  in  slope  after  70 
p Hz  at  Charlevoix  and  Llanrwst.  Also,  there  is  a  geographical  relat¬ 
ionship  evident  in  the  spectral  characteristics.  Those  measurenents 
which  are  most  strongly  influenced  by  marine  loading  effects 
(Queensbury,  Llanrwst  and  Charlevoix)  possess  the  highest  background 
noise  persisting  over  the  entire  spectrum.  It  is  apparent  that  in 
coastal  locations  the  lower  Limits  to  ground  noise  are  determined  by 
dynamical  effects  in  the  oceans.  It  is  not  clear  whether  this  is  also 
the  case  for  inland  sites. 

5.3  Comparison  of  observed  and  theoretical  tidal  tilts 

Figure  12  shows  polar  diagrams  of  the  observed  and  computed  M2  and 
Op  mean  admittances  for:  tiltmeters  105  and  107  (borehole  1);  tiltmeter 
106  (borehole  2) ;  and  the  AN\C  measurements  A/D  and  C/D  [Peters  et  al, 
1983]  made  in  the  shallow  vailt  (figure  2).  Also  shewn  are  the  loading 
model  results  based  on  the  c xnputations  outlined  in  section  4.  Peters 
et  al  [1983]  concluded  that  3ue  to  the  geometry  of  the  vault  and  the 
expected  strain-tilt  coupling  effects,  the  ANAC  measurements  are  best 
approximated  by  the  average  if  the  two  measurements.  An  error  of  180 
degrees  in  the  phase  of  the  ®iNAC  D  component  quoted  in  Peters  et  al 
[1983]  has  been  corrected. 

There  are  substantial  differences  in  the  observed  M2  results  airong 
all  of  the  measurenents.  A  3-8%  amplitude  discrepancy  in  both  direct- 


The  main  difference  between  the  spectra  is  in  the  power  associated 
with  the  4th  and  6th  diurnal  bands.  In  both  cases  the  level  in  the 
south  component  is  ten  times  smaller  than  in  the  east  component.  This 
difference  is  surprising  since  the  large  marine  tide  and  the  substant¬ 
ial  drying  areas  in  the  estuary  to  the  south  of  the  site  are  expected 
to  be  a  rich  source  of  higher  harmonics  in  the  marine  tides.  Peters  and 
Kumpel  [in  prep]  show  that  the  non-linear  tilt  signal  can  be  explained 
by  loading  of  the  marine  tides  and  that  the  large  difference  between 
the  component  directions  is  due  to  the  spatial  distribution  of  these 
shallow  water  constituents. 

It  is  interesting  to  compare  the  residual  background  noise,  after 
ranoval  of  the  tides,  among  different  measurements  which  cover  a  range 
of  techniques  and  geographical  locations.  In  figure  11,  we  compare  the 
published  background  noise  spectra  from  three  strain  and  four  tilt 
measurements,  including  that  of  107  east  frcm  borehole  1  (see  figure 
caption  for  details) .  Hie  spectra  are  plotted  on  a  log-log  scale  and 
all  show  the  characteristic  "red  noise"  linear  trend  typical  of  geoph¬ 
ysical  spectra  plotted  in  this  way.  Interestingly,  the  slopes  among  the 
strain  spectra,  and  among  three  of  the  tilt  spectra,  are  remarkably 
consistent.  In  the  case  of  strain,  the  relationship  between  noise 
power,  P,  and  frequercy,  f,  is  P<*f-2.9.  For  each  of  the  tilt  measure¬ 
ments,  P~<f~2.3,  excepting  Llanrwst  where  the  noise  level  decays  more 
slowly.  There  is  evidence  that  these  slopes  change  between  the  frequ¬ 
encies  50  and  500  ^Hz.  This  is  clear  for  the  Queensbury  strain  spect¬ 
rum,  and  also  for  the  Pinon  Flat  and  Poorman  strain  spectra  plotted  to 


eted.  The  origin  of  this  component  of  the  long  term  tilt  is  not  known. 

A  more  detailed  study  of  these  and  shorter  period  hydrologically  induc¬ 
ed  tilts  is  in  progress  [Kumpel  and  Peters,  in  prep]  using  continuous 
water  level  data.  Theoretical  results  [Kumpel,  in  prep]  suggest  that 
tilts  of  the  observed  size  can  be  induced  by  pore  pressure  changes  that 
occur  in  the  rock  matrix. 

5.2  Power  spectrum  of  the  residual  time  series 

The  accuracy  of  the  tidal  estimates  depends  upon  the  residual 
energy  in  the  tidal  bands.  Before  discussing  the  tidal  results,  we  will 
examine  the  residual  power  spectra  of  the  ensemble  average  (section 
3.3)  shown  in  figure  10  for  both  components  of  tiltmeters  106  and  107. 

The  residuals  in  all  components  share  a  number  of  common  features. 
Firstly,  there  is  a  strong  presence  of  quarter,  sixth  and  tenth  diurnal 
energy,  presumably  generated  by  shallow  water  non-linear  interactions 
from  the  marine  load  signal.  Otherwise,  the  spectra  consist  of  a 
statistically  featureless  background  which  decreases  in  power  with 
increasing  frequency  but  flattens  out  beyond  70  ;jHz  (6  cycles/day)  at  a 
level  of  0.025  nrad2/^iHz  (power  spectral  density) .  Lastly,  there  is 
a  distinct  "hole"  in  the  semi-diurnal  band,  as  expected,  where  most  of 
the  energy  has  been  removed  in  the  least  squares  analysis  for  the  tidal 
constituents.  However,  there  is  also  significant  residual  energy,  more 
prominent  in  106,  remaining  in  that  band,  suggesting  that  the  linear 
tidal  model  cannot  account  for  all  of  the  signed. 
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100  degree  azimuth  during  the  second  period  correlates  closely  with  the 
water  level  changes.  If  it  is  assuned  that  water  level  changes  cause 
the  secular  tilt,  a  40  degree  difference  in  the  direction  of  the  water 
level  effect  in  two  boreholes  80m  apart  is  not  really  surprising.  On 
the  evidence  from  drilling  operations,  the  formation  permeability  is 
likely  to  be  fracture  dominated  and  there  is  no  reason  to  expect  a 
uniform  fracture  distribution  or  orientation  over  distances  of  the 
order  of  80m. 

The  amplitude  of  the  tilt  induced  by  the  water  level  changes  is 
different  for  each  borehole.  For  105  in  borehole  1,  the  large  spring 
water  level  transient  yielded  0.4  prad  per  metre  of  water  level  change. 
Tiltmeter  107  in  the  same  borehole  gave  0.35  jurad  per  metre  for  the 
April  event  and  0.5  jurad  per  metre  for  the  November  event.  Tiltmeter 
106  in  borehole  2  was,  however,  nearly  twice  as  sensitive  on  both 
occasions,  giving  0.65  /irad  per  metre  and  0.8  prad  per  metre  in  April 
and  November,  respectively.  Again,  it  is  not  surprising  that  in  an  area 
of  fracture  dominated  permeability,  water  level  induced  tilts  are 
varying  on  a  scale  of  80m  since  the  fractures  in  the  inmediate  vicinity 
of  the  boreholes  are  expected  to  dominate  the  response. 

There  is  a  very  clear  monotonic  drift  in  107  of  1.2  /irad  to  the 
north-east  (#=  45  degrees)  in  addition  to  the  first  order  ( 0 =  135 
degrees)  cyclical  effect  that  correlates  with  the  water  level.  There  is 
evidence  of  a  similar  but  smaller  trend  in  106,  although  the  direction 
may  be  reversed.  The  sense  of  the  residual  drift  in  the  case  of  105  is 
not  certain  because  the  unnial  water  level  related  cycle  was  not  ccmpl- 
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5.  RESULTS 


5.1  Correlation  of  secular  tilt  with  water  level  changes 

Figures  8(a)  and  8(b)  show  the  low  pass  filtered  observations 
(sampled  every  48  hours)  for  tiltmeter  105  operating  in  borehole  1 
during  1980/1  and  tiltmeters  106  and  107  operating  concurrently  in 
boreholes  2  and  1  during  1981/2.  Also  shown  are  hand  measured  water 
level  data  observed  in  well  B(70)  near  borehole  2.  The  east  and  south 
components  of  105  behave  very  similarly,  both  in  the  short  and  the  long 
term.  The  total  range  of  the  105  drift  in  each  direction  over  the  260 
day  interval  is  1.5  prad.  Although  the  water  level  is  not  adequately 
sampled  to  show  short  term  fluctuations,  there  is  a  strong  long  term 
correlation  with  the  tilt.  A  similar  correlation  exists  between  water 
level  and  tilt  for  tiltmeters  106  and  107.  Furthermore,  the  direct¬ 
ions  of  highest  correlation  with  water  level  change  are  similar  for 
each  of  the  boreholes.  This  can  best  be  seen  in  the  tilt  trajectories 
plotted  in  figure  9.  The  predominant  drift  direction, 6,  for  105  has  an 
azimuth  of  135  degrees  and  for  107,  operating  in  the  same  borehole  one 
year  later,  the  drift  azimuth  is  145  degrees.  The  situation  for  106  is 
more  complicated .  We  believe  that  fran  the  time  of  installation  in  1981 
to  day  180  (June  29)  of  1982,  tiltmeter  106  was  not  well  seated  in  the 
borehole.  The  drift  signal  during  this  time  was  large  and,  apart  from 
the  large  water  level  related  transient  in  April,  bears  little  resembl¬ 
ance  to  the  107  signal.  After  day  180,  106  was  reinstalled  and  frcm 
this  time  the  drift  signal  was  considerably  smaller.  Hie  tilt  in  the 
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4.  THE  MARINE  TIDAL  LOADING  MODEL 

The  load  tilt  calculations  were  made  by  the  same  method  as  those 
of  Beaumont  and  Lambert  [1972] ,  Beaumont  [1978b] ,  and  Beaumont  and 
Boutilier  [1978] ,  assuming  a  uniform,  spherical  earth.  Firstly,  the 
Green's  functions  for  the  point  load  response  of  the  earth  model  were 
computed.  The  integrated  effect  of  the  ocean  tide  distribution  was 
found  by  convolving  the  Green's  functions  with  the  in-phase  and  quadra¬ 
ture  components  of  the  tidal  distribution.  The  cotidal  charts  were 
divided  into  triangular  areas  in  which  the  amplitude  and  phase  of  the 
tides  could  be  regarded  as  constant  in  order  that  the  convolution  could 
be  evaluated  numerically  [Bower,  1971] .  A  more  detailed  representation 
of  the  tides  in  the  St.  Lawrence  estuary  and  Saguenay  river  (figure  5) 
was  used  here  than  in  the  earlier  work  because  most  of  the  load  tilt  is 
generated  close  to  the  Charlevoix  Observatory.  The  cotidal  charts 
(figures  6  and  7)  for  this  region  are  enpirical  but  are  based  on  obser¬ 
vations  from  the  coastal  sites  indicated,  and  on  numerical  calculations 
of  the  M2  and  Oi  constituents  [El-Sabh  et  al,  1979] .  We  suspect  that 
the  largest  errors  in  the  marine  tidal  loading  model  arise  frcm  areas 
that  "dry"  during  the  tidal  cycle  (outlined  by  the  fine  dashed  lines  in 
figure  6) .  For  this  reason  the  computation  has  been  done  for  two 
extreme  cases  in  which  there  is  no  drying  and  in  which  areas  that  dry 
remain  permanently  dry.  Errors  also  arise  frcm  the  norrstationarity  and 
non-linearity  of  the  tides. 
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3.3  Fourier  Analysis  of  the  Tidal  Residuals 

Power  spectral  analysis  of  the  106  and  107  tidal  residuals  was 
performed  in  order  to  determine  an  upper  limit  of  the  signal  to  noise 
ratio  for  the  main  tidal  constituents.  The  Fourier  transform  with  the 
Hanning  window  applied  was  calculated  for  5  independent  2  month  sets 
covering  the  13  months  of  data  studied.  The  power  spectrum  was  computed 
for  each  set  and  the  ensemble  average  calculated  to  increase  the  stab¬ 
ility  of  the  estimates. 

With  a  view  to  examining  the  detailed  spectral  content  of  the 
tidal  residuals,  a  high  resolution  Fourier  amplitude  spectrum  was 
calculated,  with  the  Hanning  window  applied,  for  the  106  and  107  resid¬ 
ual  time  series,  after  subtracting  the  mean  tidal  estimates  resulting 
from  harmonic  analysis  of  the  29  overlapping  sample  series.  The  amplit¬ 
ude  spectrum  was  computed  with  a  frequency  resolution  of  58.5nHz, 
corresponding  to  the  fundamental  frequency  for  the  time  series,  at 
intervals  of  30  nHz.  Our  main  objective  in  this  high  resolution 
analysis  was  to  identify  unresolved  and  unaccounted  for  waves  which 
possess  enough  energy  to  influence  the  major  constituents  in  our  two 
months  analyses.  We  might  expect  that  the  culprits,  if  any,  will  be 
non-linear  shallow  water  waves,  for  example  OP2  or  MKS2,  which  certain¬ 
ly  do  not  have  the  same  admittance  as  their  neighbours.  In  addition, 
there  is  the  possibility  that  unresolved  astronomical  waves,  forced  to 
have  the  same  admittance  as  their  neighbours  in  a  group  under  the 
constant  admittance  assumption,  are  modelled  poorly  and  can  therefore 
modulate  one  of  the  principal  constituents. 
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of  the  results  must  take  this  into  account,  since  the  importance  of 
interaction  of  the  Hanning  data  window  with  .he  synthesised  data 
depends  on  the  position  of  the  filled  gap  wi.hin  the  sanple  epoch. 

Figure  4  shows  the  relative  weighting  of  data  by  the  Hanning  window  as 
a  function  of  position  in  the  data  window.  Data  near  the  ends  cont¬ 
ribute  little  to  the  analysis  results  when  weighted  by  the  Hanning 
window.  The  relative  contribution  of  the  synthesized  data,  in  terms  of 
its  weighting,  to  the  analysis  of  each  overlapping  sample  epoch  has 
been  computed  on  this  basis  and  each  one  classified  according  to  wheth¬ 
er  the  contribution  is  less  than  10%,  between  10%  and  20%  or  greater 
than  20%.  We  consider  that  estimates  affected  at  a  level  greater  than 
20%  are  not  reliable.  Although  this  treatment  is  not  an  ideal  solut¬ 
ion,  it  goes  further  than  most  in  dealing  with  the  difficult  problem  of 
data  continuity,  and  represents  an  attempt  to  get  the  most  out  of  what 
amounts  to  an  average  working  data  set.  More  recent  data  have  fewer 
gaps  and  their  analysis  will  not  be  contaminated  to  the  same  degree. 

A  representative  estimate  of  the  mean  tidal  admittance  is  necess¬ 
ary  for  comparison  of  the  average  coherency  among  the  three  borehole 
tiltmeters  and  the  vault  measurements,  and  to  compere  the  observations 
with  the  theoretical  tilt  based  on  a  marine  loading  model  (section  4) . 
Schuller  [1977]  argues  that  the  best,  unbiased  estimate  of  the  mean 
admittance  is  derived  f ram  the  means  of  the  n-phase  and  quadrature 
components  of  the  time  varying  admittance  function  determined  as  above, 
since  any  bias  due  to  aliasing  is  avoided.  'Che  error  in  the  estimate  is 
derived  from  the  variance  of  the  time  varyinq  admittance  function. 


equi-spaced  time  series  and  care  must  be  taken  to  avoid  aliasing.  The 
Nyquist  frequency  is  given  by  Q^Tt/at,  where  at  is  the  time  shift  of 
the  sequential  epochs.  The  maximum  frequency  present  in  the  admittance 
function  is  determined  by  the  banctoidth  of  the  Hanning  window, 

Bh=4  Tf/T,  where  T  is  the  length  of  the  saitple  epoch.  Thus  the  condit¬ 
ion,  ATCT/4,  determines  the  minimum  time  shift  required  to  avoid  alias¬ 
ing  in  the  sequential  analysis.  The  sanple  epoch  for  the  present 
analysis  was  chosen  to  be  1440  hours  in  order  to  separate  the  main 
tidal  constituents,  and  the  time  shift  set  to  288  hours,  which  satis¬ 
fies  the  above  condition.  However,  the  amplitude  of  the  first  side-lobe 
of  the  Hanning  window  is  2.6%  of  that  of  the  main  lobe,  which  therefore 
potentially  allows  sane  leakage  between  waves  separated  by  frequencies 
in  the  range  of  138.89  uHz  (4TJ/T)  to  208.33  uHz  (6  Tf/T) .  Among  the 
major  constituents,  only  N2  will  be  affected  by  significant  leakage 
(from  M2)  so  that  the  N2  admittance  function  should  be  interpreted 
carefully  on  account  of  small  but  possibly  significant  aliasing 
effects. 

Apart  from  the  fact  that  a  continuous  record  is  desirable  if  we 
hope  to  reliably  detect  departures  from  a  stable  admittance  function, 
it  is  also  logistically  necessary  if  we  wish  to  perform  the  type  of 
analysis  described  above.  Inevitably,  gaps  exist  in  our  data  and, 
after  filtering,  with  attendant  loss  of  data,  amount  to  approximately 
18%  for  106  and  12%  for  107  of  the  total  recording  period.  Gaps  can  be 
interpolated  within  the  HYCDN  procedure  using  the  synthesised  tide 
based  on  the  mean  admittance  of  the  total  data  set.  Any  interpretation 


spectral  side-lobes  associated  with  the  rectangular  window  have  an 
effect  even  in  regression  analysis.  We  are  penalised  in  this,  however, 
because  the  Hanning  spectral  window  broadens  the  spectral  peaks  by 
comparison  with  those  given  by  the  equivalent  rectangular  data  window. 
Therefore  the  minimum  data  window  must  be  increased  to  two  months  in 
order  to  avoid  the  leakage  of  energy  between  major  constituents,  espec¬ 
ially  M2  and  N2-  This  increase  correspondingly  reduces  the  number  of 
independent  estimates  of  the  admittance  for  the  second  and  longer  data 
set  to  six. 

Our  approach  to  the  analysis  of  potentially  non-stationary  tidal 
signals  is  to  assume  that  admittance  changes  are  sufficiently  slow  that 
two  month  samples  may  be  considered  to  be  stationary.  By  analysing 
samples  that  overlap  in  time  we  expect  to  detect  variations  in  the 
Earth's  tidal  response,  or  the  tidal  forcing,  as  changes  in  the  average 
admittance  from  one  sample  to  the  next.  In  practice,  the  sequence  of 
average  acknittances  (which  we  will  refer  to  as  the  admittance  function) 
may  reflect  apparent,  but  unreal  time  variations  which  result  frcm 
either  the  inaccuracy  of  the  assumption  that  the  admittance  is  constant 
within  a  wave  group,  or  from  the  effects  of  non-linear  or  other  waves 
of  a  non-astroncmical  origin  that  are  not  included  in  the  analysis.  In 
both  cases  periodicities  will  be  induced  which  result  in  a  time  varying 
acknittance  function  due  to  intermodulation.  Considerable  care  is 
required  to  properly  detect  such  modulation  and  to  distinguish  it  frcm 
real  time  variations  in  the  earth  response  or  the  tidal  forcing. 

As  pointed  out  by  Schuller  [1977] ,  the  admittance  functions  are 
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variations  in  the  constituents  where  it  is  necessary  to  avoid  the 
confusion  between  actual  changes  in  the  earth  response  or  the  tidal 
forcing,  and  apparent  admittance  changes  due  to  modulation  when  finite 
samples  of  the  data  are  analysed.  We  have  used  the  HYCDN  method  of 
Schuller  [1977]  in  which  the  original  505  astronomical  waves  have  been 
reduced  to  73  whose  theoretical  earth  tide  amplitudes  (ignoring  tidal 
loading)  are  greater  than  0.1  nr ad.  This  represents  a  compromise  bet¬ 
ween  the  cost  of  computation  and  the  precision  required  of  the 
analysis.  It  will  be  shown  Later  that  the  background  noise  is  suffic¬ 
iently  high  that  the  effects  of  omitted  astronomical  waves  are  not 
significant.  Since  waves  closely  spaced  in  frequency  cannot  be  resolved 
for  the  data  window  used,  waves  clustered  around  each  of  the  10  main 
constituents  are  grouped  anc  each  group  is  assumed  to  have  a  constant 
admittance . 

The  HYCDN  procedure  involves  the  computation  of  the  theoretical 
earth  tide,  and  the  setting  up  and  solution  of  the  least  squares  normal 
equations  in  terms  of  the  admittance  of  the  observed  data  to  the  body 
tide  input.  The  associated  error  for  the  estimates  is  computed  for  each 
tidal  band  (diurnal,  semi-d  urnal  and  ter-diumal)  on  the  basis  of  the 
frequency  dependant  residua  variance  determined  from  the  Fourier 
analysis  of  the  residual  time  series  [Schuller,  1977]. 

One  of  the  chief  advantages  of  the  HYCDN  method  is  that  it  comb¬ 
ines  the  spectral  analysis  technique  of  windowing  the  data,  using  the 
Hanning  window,  with  the  least  squares  analysis  for  the  admittance  of 
the  tidal  groups.  The  approach  recognises  that  the  effects  of  the 
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j  For  calibration  cycles  in  which  the  background  noise  was  normal,  the 

standard  deviations  we.  a  typically  at  the  level  of  the  resolution  of 
the  digitising  table  (0.1mm  or  0.5nrad).  Figure  3  shows  the  calibrat- 
I  ion  amplitudes  recorded  on  the  primary  recorders  for  tiltmeters  106  and 

107.  The  scale  factor  for  each  channel  is  determined  to  better  than 

0.1%. 

I  Prior  to  analysis  the  data  were  digitally  filtered  to  separate  the 

tides  from  the  long  period  and  secular  components  of  the  signal.  For 
the  harmonic  analysis  a  61  point  symmetrical  bandpass  filter  designed 
to  reject  the  diurnal  and  semi-diurnal  bands  was  used.  Hie  response  of 
the  filter  was  computed  for  each  tidal  frequency  in  the  harmonic 
analysis  and  the  computed  amplitudes  correct  2d  for  the  small  effects, 
f  typically  0.1%,  of  the  filtering.  A  lew  pas  ;  filter  was  used  to  pre¬ 

whiten  the  data  prior  to  the  power  and  high  resolution  spectral 
analyses  and  the  estimates  were  subsequently  corrected  for  its  effect. 

I 

3.2  Harmonic  Analysis 

The  aim  of  harmonic  analysis  of  the  til:  data  is  to  determine 
accurately  the  amplitudes  and  phases  of  cons  -ituents  of  known  frequency 
within  the  data.  For  this  study  the  constituents  of  interest  are  those 
tidal  constituents  of  astronomical  origin  which  contain  sufficient 
energy  either  to  be  recognisable  above  the  noise,  or  those  that  are 
unresolved  in  the  analysis  yet  must  be  included  in  the  solution  to 
avoid  significant  modulation  of  the  adjacent  waves  in  the  spectrum. 

The  latter  consideration  is  particularly  inportant  in  the  study  of  time 


instruments,  but  essential  for  the  study  of  time  variations  in  tidal 
admittance  from  which  the  effects  of  sensitivity  changes  of  the  tilt- 
meter  and  the  recording  system  must  be  eliminated. 

3.  DATA  ANALYSIS 

3.1  Data  Preparation 

The  data  analysed  in  this  study  were  recorded  during  the  periods 
27/Sept/8 0  to  29/June/81  (tiltmeter  105  in  borehole  1)  and  27 /Nov/81  to 
31/Dec/82  (tiltmeters  106  and  107  in  boreholes  2  and  1,  respectively). 
All  data  were  recorded  on  strip  chart  recorders  at  a  recording  speed  of 
30mm  per  hour  and  at  a  sensitvity  of  approximately  0.2  mn  per  nrad. 
Hourly  time  marks  were  superimposed  on  the  recordings  and  the  data 
digitised  at  these  marks.  Spurious  points,  resets  and  short  thermal 
trends  were  removed  and  a  few  gaps  of  under  24  hours  duration  were 
manually  interpolated  on  the  basis  of  data  before  and  after  the  gap. 
This  was  possible  because  the  long  period  trends  in  the  data  change 
slowly  in  comparison  with  the  tide. 

The  data  were  scaled  using  calibration  factors  derived  from  twice 
weekly  in  situ  calibration  sequences.  Hie  mean  and  standard  deviation 
of  each  sequence  was  computed  and,  since  most  changes  in  the  calibrat¬ 
ion  amplitudes  corresponded  with  recorder  changes,  a  scale  factor  was 
determined  for  each  recorder  channel  using  the  weighted  mean  of  the 
calibration  amplitudes  determined  from  data  registered  on  that  channel. 


lack  of  it,  in  the  tidal  constituents.  By  "true"  variability  we  mean 
fluctuations  that  require  an  interpretation  in  terms  of  time  changes  in 
the  Earth's  response  to  the  tidal  forcing,  or  changes  in  the  tidal 
forcing  itse1^  Our  aim  is  first  to  eliminate  apparent  variations  which 
are  either  artifacts  of  the  analysis  procedure,  or  are  due  to  inter¬ 
modulation  between  constituents  not  separable  over  the  analysis  epoch. 

An  examination  of  the  physical  causes  of  the  variability  that  we  estab¬ 
lish  will  be  the  subject  of  a  separate  study. 

The  results  of  the  ovei lapping  two-monthly  HYCDN  analyses  describ¬ 
ed  in  section  3.2  are  plotted  in  figures  13  to  18.  The  amplitude  and 
phase  estimates  for  tiltmeter  105  operating  in  borehole  1  for  M2,  N2, 

s2K2,  Oi  and  P.1S1K1  are  ploi  ted  against  the  central  time  of  the  data 
window  in  figures  13  and  14.  The  relative  influence  of  interpolated 
gaps  (see  section  3.2)  is  indicated  at  the  top  of  each  figure  for  each 
sample  series  by  symbols  explained  in  the  caption  for  figure  13.  The 
significance  of  admittance  variations  plotted  in  the  panels  below 
should  be  judged  in  the  light  of  these  qualitative  confidence  estimates 
because  the  effects  of  a  poor  interpolation  model  cannot  be  disting¬ 
uished  from  an  admittance  variation.  Less  weight  should  be  given  to 
estimates  from  series  for  which  the  interpolation  effect  is  great¬ 
er  than  20%  than  those  where  the  corresponding  effect  is  less  than  10%. 

In  general,  the  variab  lity  of  the  estimates  for  tiltmeter  105  is 
within  or  close  to  the  95%  confidence  envelope  (shown  by  the  dotted 
curves  in  the  plots),  indie  iting  that  fran  a  statistical  point  of  via* 
with  respect  to  the  residua  noise,  the  variations  are  hardly  signific- 
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ant.  Despite  this,  there  does  seen  to  be  an  underlying  correlation 
between  the  admittance  curves  in  the  two  component  directions  for  each 
wave.  The  M2  south  and  east  amplitudes  both  show  a  broad  peak  spanning 
the  winter  months  (December  to  March)  while  :he  phase  lags  are  reason¬ 
ably  stable.  The  south  and  east  N2  estimates  exhibit  troughs  both  in 
amplitude  and  phase  lag  during  the  winter,  while  S2  trends  steadily  to 
increasing  amplitude  and  phase  lag  through  most  of  the  9  month  period. 
Although  the  diurnal  constituents  are  not  as  well  determined,  the 
amplitude  and  phase  lag  variations  of  P]S]Ki  and  Oi,  respectively,  are 
quite  similar  in  both  azimuths. 

It  is  apparent  from  these  data  that  the  range  of  the  amplitude 
variations  is  much  the  same  for  all  of  the  constituents  (between  2  and 
6  nr ad) .  Thus,  the  fractional  amplitude  changes  in  M2  amount  to  about 
6%  and  for  the  diurnals,  up  to  50%.  We  will  look  now  at  the  results 
from  the  second  data  set  (tiltmeters  106  and  107,  figures  15  to  18)  to 
see  if  the  behaviour  noted  for  105  persists  lieyond  the  first  period  in 
borehole  1,  and  if  that  behaviour  is  coherent  between  the  two  bore¬ 
holes. 

Although  the  gaps  in  the  observed  106  aid  107  data  are  in  most 
cases  quite  short  (2  to  4  days)  and  the  number  of  gaps  small  (7  for  106 
and  8  for  107)  their  effect  has  propogated  widely  through  the  sequent¬ 
ial  analyses.  Ihe  distribution  of  gaps  is  qiite  different  for  the  two 
installations,  but  nearly  the  same  for  both  components  of  the  same 
tiltmeter.  Out  of  29  points,  only  12  are  affected  at  less  than  10%  for 
each  instrument,  26  at  less  than  20%  for  107  and  15  under  20%  for  106. 
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The  plots  therefore  have  to  be  interpreted  carefully.  Ignoring  the  gap 
problem  for  the  moment,  we  will  look  at  the  general  behaviour  of  the 
data  and  subsequently  determine  if  the  features  remain  under  closer 
scrutiny. 

In  a  similar  way  to  105  operating  in  borehole  1  the  variations  in 
the  four  components  are,  in  most  cases,  hardly  significant  in  terms  of 
the  95%  confidence  error  envelope.  Nevertheless,  among  the  semi¬ 
diurnal  waves  there  are  particularly  strong  correlations  between  the 
admittance  curves  from  the  orthogonal  components  of  the  same  instru¬ 
ments.  With  the  exception  of  N2  in  the  south,  we  can  also  detect  a 
clear  but  somewhat  weaker  correlation  between  the  estimates  determined 
in  the  same  azimuth  from  the  two  instruments.  So  it  appears,  for  the 
semi-diurnals  at  least,  that  there  is  better  agreement  between  the 
south  and  east  estimates  in  a  given  borehole  than  between  redundant 
measurements  made  in  the  two  boreholes.  This  is  disturbing  because  it 
suggests  that  a  measurement  effect  is  influencing  the  estimates.  Ihe 
patterns  of  the  diurnal  admittance  variations  (figures  17  and  18)  are, 
on  the  other  hand,  quite  distinct  in  either  direction  and  correlate 
well  between  the  boreholes. 

There  are  two  other  general  features  of  the  data  that  are  consist¬ 
ent  with  the  105  data  set.  Firstly,  the  absolute  range  of  the  variat¬ 
ions  is  similar  for  the  five  constituents  irrespective  of  wave  amplit¬ 
ude.  Secondly,  there  is  no  discernible  correlation  between  the  variat¬ 
ions  of  the  different  waves.  Thus,  we  cannot  identify  a  systenatic 
semi-diurnal  or  diurnal  characteristic  behaviour. 
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If  we  examine  the  semi-diurnal  admittance  variations  in  conjunct¬ 
ion  with  the  interpolation  distribution,  it  becomes  clear  that  for  the 
"reliable"  data  (the  subset  of  the  overlapping  analysis  results  in 
which  points  contaminated  by  interpolation  a:  a  level  greater  than  20% 
are  ignored) ,  there  is  a  strong  correlation  Ixath  between  orthogonal 
directions  of  the  same  instrument  and  fron  one  instrument  (or  borehole) 
to  another.  Moreover,  the  lack  of  correlation  between  components  from 
the  two  boreholes  can  be  explained  by  the  d: fferenoe  in  interpolation 
effects.  We  can  expect  that  if  the  south  ard  east  components  of  the 
measuranent  are  really  varying  in  the  same  way,  then  the  admittance 
curves  will  be  distorted  in  the  same  way  by  the  same  distribution  of 
interpolated  data.  Such  an  effect  would  lea/e  the  orthogonal  components 
of  each  of  the  tiltmeters  highly  correlated,  as  observed.  Since  the 
gaps  in  the  two  tiltmeter  records  occur  at  different  times,  distortions 
in  the  admittance  curves  will  arise  at  different  times  and  the  maximum 
differences  will  occur  where  the  interpolation  effect  is  large  for  one 
data  set  and  small  for  the  other.  The  major  contrasts  occur  at  analyses 
1,  9,  22,  23,  27  and  28.  Eliminating  22  and  23  of  107  M2  makes  all  four 
components  agree  in  both  amplitude  and  phase  lag.  An  improvement  is 
made  for  N2  by  disregarding  analyses  9,  22  and  23,  although  the  south 
direction  variations  remain  poorly  correlated  between  the  instruments. 

In  the  case  of  S2,  for  which  the  inter- instrument  correlations  are 
good,  little  improvement  is  gained  by  considering  interpolation 
effects. 


The  diurnal  results  are  quite  well  correlated  between  instruments, 
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yet  the  style  of  the  variations  is  different  in  the  two  azimuths.  In 
most  cases,  particularly  in  the  south  components,  the  variability  is 
greatest  during  the  winter  months  at  the  same  time  that  the  least 
squares  errors  are  largest. 

5.5  Polar  plots  of  the  variability  of  the  M2  and  Oi  admittances. 

The  results  shown  in  figures  15  to  18  can  also  be  presented  as 
polar  plots  in  which  variations  in  admittance  trace  out  a  path,  with  a 
surrounding  swath  of  error  estimates,  in  polar  coordinates.  Figures  19 
and  20  show  these  phasor  trajectory  plots  of  M2  and  Oi  for  106  and  1(77 
south  and  east  components.  The  points  are  plotted  using  symbols,  as 
described  earlier,  indicating  the  effect  of  interpolation.  Ninety-five 
per  cent  error  circles  have  been  plotted  for  each  point  and  then,  to 
avoid  confusion,  largely  erased  leaving  only  a  representative  arc 
labelled  by  the  analysis  nimber.  The  reader  can  reconstruct  the  compl¬ 
ete  plot  from  the  information  shown. 

Following  the  procedure  used  earlier  to  represent  the  weighting 
of  the  data  points  accordinq  to  the  symbols,  we  have  attempted  to 
identify  the  most  likely  pa  :hs  of  the  admittance  variations.  The  inter¬ 
pretation  is  shown  in  figure  21  which  should  be  considered  in  conjunct¬ 
ion  with  figures  19  and  20.  Looking  at  M2  first,  the  resemblance 
between  the  two  east  components  is  renarkable  and  is  firmly  based  on 
the  "under  10«"  interpolation  effect  data.  Even  the  early  data 
(analyses  1  to  5) ,  where  106  suffers  from  large  interpolation  errors, 
correlate  closely  in  their  behaviour.  The  strong  trend  fran  analyses  6 
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to  15  occurs  during  the  spring  and  comprises  a  combined  amplitude 
decrease  of  about  5%  and  phase  lag  increase  of  2.5  to  3  degrees.  From 
mid-June  to  mid-August  (analyses  15  to  20)  only  the  airplitude  changes, 
increasing  by  3.5%  in  106  and  4.5%  in  107.  This  is  followed  by  a  period 
of  greater  stability  in  analyses  21  to  25  which  are  well  determined, 
and  in  26  to  29  which  are  less  well  determined  in  tiltmeter  106. 

In  the  south  direction,  the  resemblanoe  between  the  two  M2  traj¬ 
ectories  is  less  impressive.  However  an  early,  predominantly  phase, 
trend  (analyses  1  to  5)  and  its  well  determined  reversal  during  the 
sparing  (6  to  15) ,  show  the  same  level  of  variation  as  in  the  east 
direction.  The  fluctuations  after  this  time  are  relatively  small  and 
are  not  well  correlated  between  the  instruments.  Certainly,  the  period 
covered  by  analyses  21  to  29  is  obviously  one  of  greater  stability  for 
both  components  and  it  is  probable  that  statistical  uncertainties  in 
the  analysis  results  obscure  ary  correlation  between  trajectory  paths. 

The  Op  trajectories  are  more  complex.  In  the  east  direction,  most 
of  the  variations  are  in  phase  lag  during  the  first  four  and  last  three 
months  of  the  series.  In  between,  the  admittance  is  fairly  stable.  The 
behaviour  of  the  107  south  trajectory  is  similar  to  those  of  the  east 
direction  except  that  the  variations  are  in  amplitude  rather  than  phase 
lag.  The  trajectory  of  106  south  appears  to  undergo  two  anti-clockwise 
rotations  about  the  mean  until  it  reverts  to  a  cyclical  amplitude 
variation  similar  to  107  south. 


While  these  interpretations  cannot  be  regarded  as  definitive,  they 


represent  in  our  estimation  the  most  obvious  trends  in  the  time  variat¬ 
ion  of  the  admittances.  It  is  inportant  to  note  that  by  comparison  with 
the  confidence  estimates,  only  the  maximum  excursions  are  even  margin¬ 
ally  significant  at  the  95%  level  and  that  with  only  one  instrument  we 
would  have  difficulty  in  producing  convincing  arguments  for  variabil¬ 
ity.  The  close  correlation  between  the  two  instruments  suggests  that 
our  error  estimates,  based  on  noise  estimates  over  a  whole  frequency 
band,  are  overly  pessimistic.  The  inportance  of  observations  fran  two 
or  more  instruments  is  again  confirmed. 

We  must  now  investigate  whether  the  variability  that  we  have 
described  above  is  real,  or  an  artifact  of  the  analysis  procedure. 

5.6  High  resolution  Fourier  analysis  of  the  residuals. 

The  interaction  between  the  Hanning  window  and  gaps  in  the  data 
has  already  been  discussed  and  its  effect  minimised.  The  main  concern 
now  is  the  effects  of  any  systematic  modulation  from  waves,  within  0.39 
pHz  (tidal  analysis  bandwidth)  around  the  major  constituents,  that  were 
ignored  or  poorly  modelled  n  the  harmonic  analysis. 

The  diurnal  and  semi-d  urnal  portions  of  the  high  resolution 
Fourier  spectrum  of  the  res  duals  (section  3.3)  for  both  components  of 
tiltmeters  106  and  107  are  plotted  in  figures  22  and  23  with  a  linear 
amplitude  scale.  The  "background"  noise  in  the  diurnal  band  for  all 
four  components  is  essentia. ly  flat  with  an  approximate  level  of  1 
nrad.  In  the  semi-diurnal  band,  however,  the  residual  noise  rises,  for 
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each  component,  from  a  backround  of  0.3  nr ad  to  a  broad  peak  spanning 
the  interval  21  to  23  ^iHz  with  a  particular  concentration  of  energy,  in 
the  south  components,  around  2N2.  The  broad  feature  resembles  the 
"cusp"  commonly  observed  in  marine  tides  [Munk  et  al,  1965]  and  tilt 
and  strain  observations  affected  by  marine  tidal  loading  (for  example, 
Peters  [1978];  Beavan  and  Goulty  [1977]).  It  is  attributed  (Munk  et  al, 
[1965]  to  non-linear  interactions  between  the  strong  tidal  constituents 
and  the  low  frequency  continuum,  and  shallow  water  non-linear  interact¬ 
ions  among  the  tidal  constituents  themselves.  The  cusp  rises  to  a  level 
around  0.75  nrad  in  the  east  and  1  nrad  in  the  south  and  is  consistent 
with  the  M2  constituent  FMS  errors  derived  from  the  sequential  harmon¬ 
ic  analyses  (figure  24) .  At  this  point  we  see  the  justification  for 
choosing  a  subset  of  the  original  505  astronomical  waves  for  the  HYGON 
analysis.  We  chose  a  cutoff  for  the  body  tide  amplitude  of  0.1  nrad 
(leaving  73  waves  in  the  analysis)  and  few,  if  any,  of  the  neglected 
waves  appear  to  be  significant  in  the  residual  spectrum.  These  constit¬ 
uents  are  sufficiently  below  the  observed  background  noise  to  make 
inter-constituent  modulation  insignificant,  even  for  a  signal  augmented 
by  marine  loading.  The  effect  of  the  general  background  noise  level  has 
already  been  included  in  the  error  analysis  for  the  two  monthly  data 
sets. 


A  closer  inspection  of  the  semi-diurnal  spectrum  shows  holes  in 
the  vicinity  of  the  main  tidal  lines  which  results  fran  the  least 
squares  removal  of  energy  at  those  frequencies.  Rising  significantly 
above  the  background  in  all  four  components,  but  most  pronounced  in 


106,  is  a  pair  of  peaks  located  at  a  distance  of  one  resolution  band¬ 
width  either  side  of  M2.  This  is  the  spectral  representation  of  the 
observed,  "around  6  month"  cycle  identified  in  the  M2  admittance  varia¬ 
tions.  That  this  is  an  isolated  and  significant  pair  indicates  that  the 
variation  is  real  and  not  due  to  intermodulation  effects.  No  such 
behaviour  is  apparent  for  N2,  or  S2  which  showed  a  less  significant 
variation  in  amplitude.  The  residual  spectra  of  both  azimuths  of  107 
have  a  potentially  significant  line  near  22  pHz  which  may  introduce  a 
degree  of  modulation.  The  source  of  the  60  to  80  day  periodicity  noted 
in  N2  is  not  clear,  although,  given  the  high  residual  amplitude  levels 
near  the  N2  frequency,  sane  modulation  may  be  expected. 

The  two  prominent  peaks  in  the  south  direction  (near  21.5  ^aHz) 
coincide  with  the  2N2  and  p2  constituents  which  are  included  in  the 
linear  harmonic  analysis.  Cod in  [1973],  in  his  analysis  of  marine 
tides  at  Quebec,  has,  however,  identified  p 2  as  a  double  constituent, 
sharing  its  frequency  with  the  non-linear  interaction  2M2-S2.  He  has 
also  detected  the  nonlinear  wave  O2  between  2N2  and  p2.  A  similar 
effect  may  occur  in  the  loading  near  the  Charlevoix  site.  If  this  is 
so,  the  non-linear  part  will  remain  because  the  2N2  group  cannot  be 
realistically  modelled  with  a  constant  admittance. 

6.  DISCUSSION 

Spatial  coherency  of  tit  observations  is  essential  if  they  are  to 
be  interpreted  in  terms  of  regional  crustal  properties  and  processes. 
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The  general  lack  of  spatial  coherency  among  measurements  within  a 
region  (or  even  beside  one  another)  is  related  to  the  ubiquitous 
presence  of  meteorological  noise,  and  the  departure  on  all  length 
scales  of  the  crust  from  lateral  homogeneity.  Ihe  advent  of  borehole 
and  long  baseline  measurements  has  increased  the  potential  for  achiev¬ 
ing  interpretable  results.  However,  most  of  the  traditional  objectives 
of  earth  tide  research  -  refining  global  Love  numbers,  discriminating 
between  crustal  models  and  solving  the  inverse  marine  loading  problem  - 
require,  as  Baker  [1979]  points  out,  regionally  representative  tidal 
response  estimates  accurate  to  better  than  2-5% .  It  is  considered 
unlikely  that  locally  generated  tilt  perturbations  and  associated  model 
corrections  will,  in  general,  be  reduced  to  those  levels.  Our  experi¬ 
ence  at  Charlevoix  provides  a  dramatic  demonstration  of  the  problem. 

The  borehole  installations  were  designed  with  a  view  to  minimising  the 
sensitivity  of  the  measurements  to  local  strain-tilt  coupling  effects. 
Yet  the  mean  M2  tidal  anplitudes  differ  by  20%  between  borehole  1  and 
borehole  2  with  no  observable  discontinuity  present  in  the  terrain. 
Effects  of  this  magnitude  are  rare  but  disturbing  and  serve  to  emphas¬ 
ise  the  difficulties  associated  with  determining  absolute  response 
estimates  with  any  confidence . 

Recent  crustal  deformation  research  using  the  earth  tide  has 
focused  on  the  investigation  of  large-scale  crustal  discontinuities 
such  as  batholiths  and  thrust  faults  [Edge  et  al,  1983;  Bonatz  et  al , 
1983] .  So  long  as  a  reasonable  degree  of  spatial  coherency  can  be 
established,  the  large  (up  to  50%)  anomalies  generated  by  lateral 


page  33 

elastic  property  contrasts  may  be  able  to  provide  constraints  on  the 
structure  of  these  features  that  are  complementary  to  seismic  studies 
[Harrison,  1976] . 

In  studies  of  time  variations  in  the  tidal  admittance  we  are  no 
longer  shackled  to  the  need  for  a  structurally  interpretable  measure¬ 
ment.  There  remains,  however,  the  need  for  spatial  coherency  of  the 
relative  changes  in  admittance  over  length  scales  smaller  than  the 
prospective  active  area.  It  is  therefore  essential,  as  in  all  regional 
studies  using  point  observations,  that  redundant  measurements  be  made. 

We  have  observed  remarkable  agreement  between  the  time  varying  admitt¬ 
ances  determined  in  two  boreholes,  lending  credence  to  a  regional 
interpretation.  Nevertheless,  the  study  of  relative  changes  in  the 
tidal  parameters  introduces  a  new  set  of  interpretation  problems  and 
ambiguities. 

There  are  few  acknowledgements  in  published  studies  of  time  vary¬ 
ing  tidal  admittance  of  the  inherent  difficulties  in  fitting  a  station¬ 
ary,  deterministic  regression  model  to  a  finite  length,  gappy  repres¬ 
entation  of  a  non-stationary,  partly  deterministic  collection  of  proc¬ 
esses.  Schuller  [1977]  has  made  an  inportant  contribution  to  the  treat¬ 
ment  of  this  problem,  allowing  for  a  "piecewise  stationary"  model  and 
using  the  Hanning  data  window  to  avoid  modulation  due  to  spectral 
leakage.  Complications  arising  from  the  interaction  of  interpolated 
gaps  with  the  data  window  have  been  addressed  in  the  present  study  as 
well  as  the  importance  of  verifying  the  absence  of  intermodulation 
through  detailed  examination  of  the  mean  admittance  residual  spectrum. 


Interpretation  of  the  determined  semi-diurnal  and  diurnal  admitt¬ 
ance  variations  at  Charlevoix  is  in  progress.  Since  the  major  component 
of  the  signal  at  the  Charlevoix  site  arises  from  loading  by  the  St. 
Lawrence  River  tide,  we  expect  that  variations  in  the  load  will  be 
reflected  in  the  tilt.  Indeed,  as  noted  in  section  5,  the  observed 
tilt  admittance  variations  neither  indicate  consistent  fractional 
changes  in  amplitude,  nor  characteristic  diurnal  or  semi-diurnal 
patterns,  both  features  that  are  expected  to  accompany  changes  in 
crustal  response.  Preliminary  estimates  of  the  M2  admittance  variations 
from  the  tide  gauge  at  St.  Joseph  de  la  Rive  immediately  south  of  the 
site  correlate  reasonably  well  with  the  tilt  results.  Thus  it  is  clear 
that  a  detailed  evaluation  of  the  stability  of  the  input  must  be  made 
before  conclusions  can  be  reached  about  possible  changes  in  crustal 
response . 

The  prominent  features  in  the  borehole  secular  tilt  correlate  with 
the  transient  and  seasonal  behaviour  of  the  water  table.  Although  the 
tilt  series  are  not  long  enough  to  indicate  established  long  term 
trends,  we  can  get  an  idea  of  the  lower  limit  of  the  detectability  of 
regional  signals  from  the  difference  in  the  long  term  drift  rates 
between  the  redundant  106  and  107  measurements.  A  crude  estimate, 
based  on  the  5  months  of  data  following  the  reinstallation  in  July, 

1982  of  tiltmeter  106,  is  2  ^rad/yr  in  the  south  direction  and  0.4 
^jrad/yr  in  the  east.  Ihe  difference  between  these  rates  points  to 
possible  problems  with  reference  point  stability.  By  comparison, 
Buchbinder  et  al  [in  press]  conclude  that,  after  the  reduction  of 
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groundwater  effects  in  the  '  0m  baseline  levelling  array  at  the 
Charlevoix  site,  the  uncertainty  with  which  tectonic  tilts  could  be 
detected  falls  withir.  a  rare  e  of  4-6  jar  ad.  In  terms  of  drift  rate, 
their  data  suggest  that  a  persistent  trend  over  a  year  or  more  of  4 
/irad/yr  would  be  detectable.  This  is  considerably  worse  than  the 
borehole  results  and  arises  frcm  the  dominance  of  groundwater  (and 
probably  thermal)  effects  which  are  at  least  5  times  larger  than  in  the 
boreholes.  During  the  comparative  tilt  experiment  at  Pinon  Flat,  Wyatt 
et  al  [1982]  measured  long  term  drift  rates  of  approximately  4  prad/yr 
in  the  26m  borehole  and  0.7  prad/yr  using  the  surface  535m  baseline 
fluid  tiltmeter.  These  results  suggest  that  the  baseline  of  a  surface 
tiltmeter  needs  to  be  an  order  of  magnituoe  greater  than  the  borehole 
depth  to  achieve  comparable  performance. 

7 .  CONCLUSIONS 


An  array  of  three  borehole  tiltmeters  has  been  established  at 
Charlevoix  Observatory,  Quebec  to  study  the  tidal  and  secular  response 
of  the  crust  within  the  Charlevoix  seismic  zone.  Differences  in  the  Mo 
tidal  results  as  large  as  20%  in  amplitude  and  5  degrees  in  phase  were 
found  between  two  boreholes  80m  apart.  Theoretical  M2  amplitudes  and 
phases  computed  from  the  tidal  loading  model  of  the  St.  Lawrenoe  River 
were  systematically  different  frcm  all  of  the  observed  results,  indic¬ 
ating  that  unconsidered  local  topographic  and  geological  strain-tilt 
coupling  effects,  in  addition  to  regional  structural  features  are  to 
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seme  extent  distorting  the  response.  The  0]  mean  admittances  were 
consistent  among  the  different  observations  and  in  good  agreement  with 
theory. 

Time  variations  in  the  tidal  admittance  were  observed  for  the 
major  constituents  M2,  N2,  S2K2,  01  and  PIS1K1.  Although  the  variat¬ 
ions  were  not  in  general  statistically  significant  at  the  95%  level 
(based  on  the  least  squares  residual  variance)  they  were  ranarkably 
consistent  between  the  two  boreholes,  indicating  that  errors  based  on 
the  residual  variance  tend  to  overestimate  the  uncertainty  of  the 
admittance  determination.  It  is  believed  that  fluctuations  in  the  load 
tide  are  at  least  partly  responsible  for  the  observed  behaviour.  How¬ 
ever,  it  is  important  to  note  that  the  physiol  significance  of  the 
variations  would  have  remained  undetected  if  there  had  been  measure¬ 
ments  only  from  a  single  borehole.  Once  again  the  importance  of  red¬ 
undant  measurements  is  confirmed. 

The  borehole  secular  tilt  is  correlated  with  water  table  fluctuat¬ 
ions,  the  typical  coefficient  being  approximately  0.5  ;urad/m.  A  prelim¬ 
inary  estimate  of  the  lower  detectable  limit  of  long  term  regional 
ananalies,  using  simultaneous  measuranents  fccm  two  boreholes,  is  0.4 
^rad/yr  -  comparable  with  the  lowest  rates  s  >  far  reported. 
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Figure  23.  Seni-diurnal  band  high  resolution  Fourier  amplitude  residual  spectra 
of  the  106  and  107  south  and  east  components.  Details  as  described  in 
figure  22. 

Figure  24.  Plots  versus  time  of  the  106  and  107  south  and  east  RMS  amplitude 
errors  associated  with  M2  amplitude  estimates  fran  the  HYODN 
sequential  analysis. 
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lags  for  0^  and  P1S1K1.  Details  as  described  for  figure  13. 

Figure  19.  Polar  trajectory  plots  of  the  106  and  107  south  and  east  M2 

admittances  estimated  f ran  the  HYCDN  sequential  analysis.  The 
observed  mean  admittance  vector  is  shewn  in  the  upper  part  of  each 
diagram.  The  shape  surrounding  the  end  point  of  the  vector 
corresponds  with  the  dashed  curve  surrounding  the  detailed  trajectory 
plot.  Each  point  on  the  trajectory  is  represented  by  a  symbol 
indicating  the  degree  to  which  that  estimate  is  influenced  by 
interpolated  gaps  (see  text) .  Filled  circles  indicate  an  influence 
less  than  10%,  squares  less  than  20%  and  triangles  greater  than  20%. 
Points  associated  with  triangles  are  least  reliable.  Part  of  the  95% 
error  circle  associated  with  each  estimate  is  indicated  by  a  matching 
sequence  number  in  sloping  numerals. 

Figure  20.  Polar  trajectory  plots  of  the  106  and  107  south  and  east  Op 

admittances  estimated  from  the  HYG>N  sequential  analysis.  Details  as 
described  in  figure  19. 

Figure  21.  Polar  trajectory  interpretation.  (Xialitative  representation  of  the 
admittance  time  variations  shown  in  figures  IS  and  20. 

Figure  22.  Diurnal  band  high  resolution  Fourier  amplitude  residual  spectra  of  the 
106  and  107  south  and  east  components,  based  on  the  mean  tidal 
admittance  determined  from  the  HYOiN  sequential  analysis.  Spectral 
estimates  are  plotted  at  approximai  ely  half  the  resolution  bandwidth 
of  58.5nHz. 
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Figure  13.  Plots  versus  time  of  the  south  and  east  tidal  amplitudes  for  M2,  N2, 
s2>  Oi  and  PlSiKi  estimated  fran  the  HYCDN  sequential  analysis  of 
tiltmeter  105  data  fran  borehole  1.  The  solid  curves  join  points 
representing  the  amplitudes  while  the  dotted  curves  indicate  the  95% 
error  envelopes.  The  points  are  positioned  in  time  at  the  middle  of 
each  analysis  epoch.  Symbols  at  the  top  of  the  diagram  indicate  the 
degree  to  which  the  individual  estimates  are  influenced  by 
interpolated  gaps  (see  text) .  Filled  circles  indicate  ar  influence 
less  than  10%,  squares  between  10%  and  20%  and  triangles  greater  than 
20%.  Points  associated  with  triangles  are  least  reliable. 

Figure  14.  Plots  versus  time  of  the  south  and  east  Greenwich  phase  lags  for 

M2fN2,S2,Oi  and  PlSiKi  estimated  fran  the  HYCDN  sequential  analysis  of 
tiltmeter  105  data  fran  borehole  1.  Details  as  described  for  figure 
13. 

Figure  15.  Plots  versus  time  of  the  south  and  east  tidal  anplitudes  for  M2,N2  and 
S2  estimated  fran  the  HYCDN  sequential  analysis  of  tiltmeters  106  and 
107  data  from  simultaneous  operation  in  boreholes  2  and  1, 
respectively.  Details  as  described  for  figure  13. 

Figure  16.  Plots  versus  time  of  the  106  and  107  south  and  east  Greenwich  phase 
lags  for  M2,N2  and  S2.  Details  as  described  for  figure  13. 

Figure  17.  Plots  versus  time  of  the  106  and  107  south  and  east  tidal  amplitudes 
for  Oi  and  PlSiKi.  Details  as  described  for  figure  13. 

Figure  18.  Plots  versus  time  of  the  106  and  107  south  and  east  Greenwich  phase 


arrows.  The  error  bars  represent  95%  confidence  bounds  on  the 
randomness  of  the  background  noise.  Logarithms  are  to  base  10. 

Figure  11.  Comparison  plot  of  tilt  and  strain  background  noise  spectra  taken  fran 
various  sites:-  Pinon  Flat,  California:  laser  strainmeter  [Berger  and 
Levine,  1974],  fluid  tiltmeter  [Wyatt  et  al,  1982];  Queensbury,  U.K.: 
laser  strainmeter  [Beavan  and  Goulty,  1977];  Poorman  mine,  Colorado: 
laser  strainmeter  [Berger  and  Levine,  1974],  Hughes  bubble  tiltmeter 
[Harrison,  1976];  Llanrwst,  U.K..  fluid  tiltmeter  [Peters,  1978] ; 
Charlevoix,  Quebec:  borehole  tiltmeter  (present  study).  The  sampled 
points  are  represented  by  filled  ci rcles  in  the  case  of  strain  and 
filled  triangles  for  the  tilt  measurements.  Logarithms  are  to  base 
10. 

Figure  12.  Polar  diagrams  showing  the  comparison  with  theory  of  the  mean  M2  and 
Oi  admittances  observed  for  the  south  and  east  components  of 
tiltmeters  105  (borehole  1),  106  (borehole  2),  107  (borehole  1)  and 
the  ANAC  tiltmeters  operating  in  the  observatory  vault.  G  is  the 
Greenwich  phase  lag.  The  ANAC  results  are  taken  from  Peters  et  al, 
[1983]  and  rotated  into  their  south  and  east  components  by  combining 
the  redundant  A  and  C  measurements  with  the  single  orthogonal  D 
measurement.  The  theory  is  represented  in  the  upper  part  of  each 
diagram  by  the  sum  of  the  body  tide  (vector  1) ,  the  far-field  load 
(vector  2)  and  the  local  load  (vectors  3  and  4;  see  text).  The  lower 
part  of  the  diagrams  shews  the  detailed  comparison  between  the 
observations  and  theory.  The  circlas  represent  95%  errors  on  the  mean 
admittance  estimates. 
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the  Charlevoix  site  shown  below  it  in  detail.  It  is  the  loading  from 
the  estuary  that  contributes  the  local  or  '3'  and  '4'  components  of 
the  tilt  phasors  in  figure  12. 

Bnpirical  M2  00 tidal  chart  for  the  St.  Lawrence  estuary,  with  detail 
of  the  area  adjacent  to  the  Charlevoix  site.  Dots  represent  the 
location  of  tide  gauge  installations. 

7.  Bnpirical  Ox  cotidal  chart  for  the  St.  Lawrence  estuary,  with  detail 
of  the  area  adjacent  to  the  Charlevoix  site.  Dots  represent  the 
location  of  tide  gauge  installations. 

8.  Plots  versus  time  of  secular  tilt  for  a)  tiltmeter  105  operating  in 
borehole  1  from  Sept.  1980  to  June,  1981;  b)  tiltmeters  106  and  107 
operating  in  boreholes  2  and  1,  respectively,  fran  Nov.,  1981  to  Dec., 
1982.  Also  shown  for  both  periods  are  manual  measurements  of  water 
level  taken  from  well  B(70) . 

9.  Secular  tilt  trajectories  of  tiltmeter  105  operating  in  borehole  1  and 
tiltmeters  106  and  107  operating  in  boreholes  2  and  1  for  the  same 
periods  indicated  in  figure  8.  Each  trajectory  describes  the  path 
taken  by  the  bottom  of  the  freely  hanging  pendulum  relative  to  the 
instrument  body.  Reference  times  are  marked  by  full  circles  and  show 
the  Julian  day  number  and  year. 

10.  Residual  power  spectra  of  the  south  and  east  components  of  tiltmeters 
106  and  107  operating  in  boreholes  2  and  1,  respectively.  Ihe 
observed  M2  and  Ox  power  levels  are  indicated  by  the  tip  of  the 
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CAPTIONS 

1.  Location  of  the  Charlevoix  Observatory.  The  lower  inset  shows  the 
observatory  area  in  relation  to  the  St.  Lawrence  Valley.  The  upper 
inset  shows  the  spatial  relationship  of  seismicity  in  the  region  to 
the  zone  of  deformation  associated  with  the  Charlevoix  crater  (after 
Lyons  et  al  [1980]). 

2.  Upper  part.  Plan  view  of  the  Charlevoix  Observatory  showing  the 
locations  of  the  three  boreholes;  water  wells  A(70  m  deep),  B(two 
wells  30  m  and  70  m  deep)  and  C(134  m  deep) ;  the  tilt  and  strain 
observatory  vault;  and  the  40  m  baseline  levelling  array.  Lower  part. 
Cross-sectional  view  of  the  borehole  layout. 

3.  Plot  versus  time  of  calibration  amplitudes  for  recorder  channels 
monitoring  tiltmeter  components  106X,  106Y,  107X  and  107Y.  Each 
amplitude  determination  is  represented  by  a  vertical  bar  the  mid-point 
of  which  is  the  amplitude,  and  the  length  equal  to  twice  the  95% 
error.  The  amplitudes  are  normalised  by  the  weighted  mean  of  each 
sequence. 

4.  Estimating  the  weighting  of  a  data  subset  as  a  function  of  its 
position  within  the  Hanning  data  window.  For  example,  a  data  subset 
spanning  the  interval  from  0.8T  to  0.9T  will  contribute  4.1%  of  the 
total  weight. 

5.  Triangular  division,  for  the  load  calculations,  of  the  whole  of  the 


St.  Lawrence  River  estuary  (upper  diagram)  with  the  area  adjacent  to 
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Figure  22 
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